SLAC-PUB-9285 
September 2002 



Search for time-dependent B® — B® oscillations 

using a vertex charge dipole technique 

Kenji Abe/ 15 ) Koya Abe/ 24 ) T. Abe/ 21 ) I. Adam/ 21 ) H. Akimoto/ 21 ) D. Aston/ 21 ) 
K.G. Baird/ 11 ) C. Baltay/ 30 ) H.R. Band/ 29 ) T.L. Barklow/ 21 ) J.M. Bauer/ 12 ) 

G. Bellodi/ 17 ) R. Berger/ 21 ) G. Blaylock/ 11 ) J.R. Bogart/ 21 ) G.R. Bower/ 21 ) J.E. Brau/ 16 ) 
M. Breidenbach/ 21 ) W.M. Bugg/ 23 ) T.H. Burnett/ 28 ) P.N. Burrows/ 17 ) A. Calcaterra/ 8 ) 

R. Cassell/ 21 ) A. Chou/ 21 ) H.O. Cohn/ 23 ) J.A. Coller/ 4 ) M.R. Convery/ 21 ) R.F. Cowan/ 13 ) 

G. Crawford/ 21 ) C.J.S. Damerell/ 19 ) M. Daoudi/ 21 ) N. de Groot/ 2 ) R. de Sangro/ 8 ) 
m '■ D.N. Dong/ 21 ) M. Doser/ 21 ) R. Dubois/ 21 ) I. Erofeeva/ 14 ) V. Eschenburg/ 12 ) S. Fahey/ 5 ) 
(N| '■ D. Falciai/ 8 ) J.R Fernandez/ 26 ) K. Flood/ 11 ) R. Frey/ 16 ) E.L. Hart/ 23 ) K. Hasuko/ 24 ) 

S.S. Hertzbach/ 11 ) M.E. Huffer/ 21 ) M. Iwasaki/ 16 ) D.J. Jackson/ 19 ) P. Jacques/ 20 ) 
J.A. Jaros/ 21 ) Z.Y. Jiang/ 21 ) A.S. Johnson/ 21 ) J.R. Johnson/ 29 ) R. Kajikawa/ 15 ) 
M. Kalelkar/ 20 ) H.J. Kang/ 20 ) R.R. Kofler/ 11 ) R.S. Kroeger/ 12 ) M. Langston/ 16 ) 
D.W.G. Leith/ 21 ) V. Lia/ 13 ) C. Lin/ 11 ) G. Mancinelli/ 20 ) S. Manly/ 30 ) G. Mantovani/ 18 ) 
T.W. Markiewicz/ 21 ) T. Maruyama/ 21 ) A.K. McKemey/ 3 ) R. Messner/ 21 ) K.C. Moffeit/ 21 ) 
T.B. Moore/ 11 ) M. Morii/ 21 ) D. Muller/ 21 ) V. Murzin/ 14 ) S. Narita/ 24 ) U. Nauenberg/ 5 ) 

H. Neal/ 30 ) G. Nesom/ 17 ) N. Oishi/ 15 ) D. Onoprienko/ 23 ) R.S. Panvini/ 27 ) C.H. Park/ 22 ) 

I. Peruzzi/ 8 ) M. Piccolo/ 8 ) L. Piemontese/ 7 ) R.J. Piano/ 20 ) R. Prepost/ 29 ) 
C.Y. Prescott/ 21 ) B.N. Ratcliff/ 21 ) J. Reidy/ 12 ) PL. Reinertsen/ 26 ) L.S. Rochester/ 21 ) 
PC. Rowson/ 21 ) J.J. Russell/ 21 ) O.H. Saxton/ 21 ) T. Schalk/ 26 ) B.A. Schumm/ 26 ) 
J. Schwiening/ 21 ) V.V. Serbo/ 21) G. Shapiro/ 10 ) N.B. Sinev/ 16 ) J.A. Snyder/ 30 ) 

H. Staengle/ 11 ) A. Stahl/ 21 ) P. Stamer/ 20 ) H. Steiner/ 10 ) D. Su/ 21 ) F. Suekane/ 24 ) 
A. Sugiyama/ 15) S. Suzuki/ 15 ) M. Swartz/ 9 ) F.E. Taylor/ 13 ) J. Thorn/ 21 ) E. Torrence/ 13) 
T. Usher/ 21 ) J. Va'vra/ 21 ) R. Verdier/ 13 ) D.L. Wagner/ 5 ) A.P Waite/ 21 ) S. Walston/ 16 ) 
A.W. Weidemann/ 23 ) E.R. Weiss/ 28) J.S. Whitaker/ 4 ) S.H. Williams/ 21 ) S. Willocq/ n) 

R.J. Wilson/ 6 ) W.J. Wisniewski/ 21 ) J.L. Wittlin/ 11 ) M. Woods/ 21 ) T.R. Wright/ 29 ) 
R.K. Yamamoto/ 13) J. Yashima/ 24) S.J. Yellin/ 25 ) C.C. Young/ 21 ) H. Yuta. (1) 



(The SLD Collaboration) 



^'Aomori University, Aomori , 030 Japan, 
( 2 ) University of Bristol, Bristol, United Kingdom, 
^Brunei University, Uxbridge, Middlesex, UB8 3PH United Kingdom, 



1 



^Boston University, Boston, Massachusetts 02215, 
^ University of Colorado, Boulder, Colorado 80309, 
^Colorado State University, Ft. Collins, Colorado 80523, 
^INFN Sezione di Ferrara and Universita di Ferrara, 1-44100 Ferrara, Italy, 
&INFN Lab. Nazionali di Frascati, 1-00044 Frascati, Italy, 
^ Johns Hopkins University, Baltimore, Maryland 21218-2686, 
^Lawrence Berkeley Laboratory, University of California, Berkeley, California 94720, 
( n ) University of Massachusetts, Amherst, Massachusetts 01003, 
( 12 ) University of Mississippi, University, Mississippi 38677, 
^Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, 
^Institute of Nuclear Physics, Moscow State University, 119899, Moscow Russia, 
( 15 ) Nagoya University, Chikusa-ku, Nagoya, 4®4 Japan, 
(16) u n i vers ity f Oregon, Eugene, Oregon 97403, 
Oxford University, Oxford, 0X1 3RH, United Kingdom, 
( 18 ^INFN Sezione di Perugia and Universita di Perugia, 1-06100 Perugia, Italy, 
^Rutherford Appleton Laboratory, Chilton, Didcot, Oxon 0X11 OQX United Kingdom, 

^Rutgers University, Piscataway, New Jersey 08855, 
( 21 ) Stanford Linear Accelerator Center, Stanford University, Stanford, California 94309, 
( 22 ) Soongsil University, Seoul, Korea 156-743, 
(23) University f Tennessee, Knoxville, Tennessee 37996, 
( 24 ) Tohoku University, Sendai 980, Japan, 
^University of California at Santa Barbara, Santa Barbara, California 93106, 
( 26 ) University of California at Santa Cruz, Santa Cruz, California 95064, 
( 2? ) Vanderbilt University, Nashville, Tennessee 37235, 
^University of Washington, Seattle, Washington 98105, 
(29) jj n i vers ity f Wisconsin, Madison, Wisconsin 53706, 
( 30 ) Yale University, New Haven, Connecticut 06511. 



(Dated: February 7, 2008) 



2 



Abstract 

We report a search for B® — B® oscillations using a sample of 400 000 hadronic Z° decays collected 
by the SLD experiment. The analysis takes advantage of the electron beam polarization as well 
as information from the hemisphere opposite that of the reconstructed B decay to tag the B 
production flavor. The excellent resolution provided by the pixel CCD vertex detector is exploited 
to cleanly reconstruct both B and cascade D decay vertices, and tag the B decay flavor from 
the charge difference between them. We exclude the following values of the B® — B® oscillation 
frequency: Am s < 4.9 ps _1 and 7.9 < Am s < 10.3 ps _1 at the 95% confidence level. 

PACS numbers: 13.20.He, 13.25.Hw, 14.40.Nd 
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I. INTRODUCTION 



Neutral B meson mixing plays a crucial role in the determination of the Cabibbo- 
Kobayashi-Maskawa matrix elements V ts and V t d- However, the extraction of these elements 
from measurements of the oscillation frequency is complicated by the presence of hadronic 

n 

uncertainties. Lattice QCD calculations [l| give a 15-20% uncertainty in the determina- 
tion of \Vtd\ from the precisely measured value of the B® oscillation frequency, Am^, due 
to the large uncertainty in the decay constant and the "B" parameter for B® mesons. This 
uncertainty can be reduced by extracting the ratio |V^ s /Vtd| from the ratio of B® and B® 
oscillation frequencies, Arris/Ama, as many theoretical uncertainties common to B° s and B® 
mixing cancel, leading to a reduced theoretical uncertainty of 5-10% . In the frame- 
work of the Standard Model, \Vt s \ is constrained by unitarity so the measurement of the B® 
oscillation frequency will significantly improve our knowledge of \Vtd\- 

Measuring the oscillation frequency requires three ingredients: the B° or B° flavor at 
both production and decay, and the proper decay time. In the Standard Model, one expects 
Am s to be an order of magnitude larger than Am<2, making it difficult to measure as the 
rapid oscillations have to be resolved in the detector. The SLD pixel CCD vertex detector 
is particularly well suited to this task thanks to its excellent three-dimensional vertex res- 
olution. The analysis presented here determines the B flavor at production by exploiting 
the large forward-backward asymmetry of polarized Z° — > bb decays and uses additional 
information from the hemisphere opposite that of the reconstructed B decay. The B flavor 
at decay is tagged by the charge difference between the B and cascade D decay vertices. 
This novel "charge dipole" technique relies heavily on the high resolution of the vertex de- 
tector to reconstruct separate secondary and tertiary vertices originating from the B — > D 
decay chain. Throughout this paper, when reference is made to a specific state, the charge 
conjugate state is also implied. 



II. APPARATUS 



The analysis uses a sample of 400 000 hadronic Z° decays collected by the SLD experiment 
at the SLC between 1996 and 1998. The detector elements most relevant to this analysis 

n 

are the pixel CCD Vertex Detector (VXD) [3( for precise track position measurements near 
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the SLC Interaction Point (IP), and the Central Drift Chamber (CDC) for charged particle 
reconstruction and momentum measurement. Charged tracks are reconstructed using hits in 
both VXD and CDC. The track impact parameter resolution at high momentum is measured 
to be 7.8 /im transverse to the beam direction (xy plane) and 9.7 fim in the plane containing 
the beam direction (rz plane). The centroid of the stable, micron-sized IP in the xy plane is 
reconstructed using tracks in sets of ~ 30 sequential hadronic Z° decays, with a measured 
precision of o/p = 3.5 /im. The z coordinate of the IP is determined event by event using 
the median z position of the tracks at their point of closest approach to the beam line. A 
precision of o z ~ 17 /im is achieved for Z° —> bb events, as estimated from the Monte Carlo 
(MC) simulation. 

For a description of the SLD detector and the MC simulation see Ref . |4j . Decays of B Q S 
mesons are modeled according to the decay modes of B® mesons [5] assuming SU(3) flavor 
symmetry. The simulation has been tuned to reproduce tracking efficiencies and impact 
parameter resolutions measured in the data. 



III. EVENT SELECTION 

The selection of 6-hadron candidates proceeds in two stages. First, hadronic Z° decays are 
selected. Second, Z° — > bb events are selected with an inclusive topological reconstruction 
of 6-hadron decays. Vertex and kinematical information (mass and momentum) is used at 
this stage to remove udsc background. This second stage also serves as a means to select 
tracks associated with the 6-hadron decay chain. 

In the first stage, we select a sample of 310 488 hadronic Z° decays. The criteria, detailed 



in Ref. aim to remove leptonic final states, two-photon collisions and beam-induced 
background events. The remaining background, predominately due to Z° — > t + t - events, 
is estimated to be < 0.1%. 

In the second stage, each event is divided into two hemispheres with respect to the thrust 
axis. In each hemisphere, a search is made for three-dimensional space points with high 
track overlap density which are displaced from the IP, taking the individual track resolution 
functions into account p. At this stage, no attempt is made to separate B and D decay 
vertices and a single "seed" vertex (SV) is formed containing all tracks from the 6-hadron 
decay chain. Tracks are required to have > 3 VXD hits, momentum transverse to the beam 
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line p± > 250 MeV/c, and three-dimensional impact parameter < 3 mm. Tracks consistent 
with originating from a 7 conversion or from K° or A decay are removed. Identified vertices 
are required to be within a radius of 2.3 cm of the center of the beam pipe to remove vertices 
resulting from interactions with the detector material. Two-prong vertices are required to 
have an invariant mass at least 0.015 GeV/c 2 away from the nominal K® mass to remove 
most of the remaining K® decays. A set of two neural networks is used to suppress light- 
flavor iudsc) background and to select the charged tracks associated with the 6-hadron 
decay |7[. The first neural network takes as input the distance between the IP and the SV, 
that distance normalized by its error, and the angle between the vertex axis and the total 
momentum vector of the SV. The vertex axis is defined as the line joining the IP and the 
SV. At least one good vertex is found in 72.7% of bottom, 28.2% of charm, and 0.41% of uds 
quark hemispheres in the MC simulation. The second neural network improves the 6-hadron 
charge and mass reconstruction by associating tracks not used in the initial vertexing to 
the seed vertex. For each track, five parameters are taken as input to the neural network. 
The first four (measured at the point of closest approach between the track and the vertex 
axis) are the distance of closest approach, the distance between the IP and the point of 
closest approach along the vertex axis, the ratio between that distance and the IP-to-SV 
distance, and the angle between the track and the vertex axis. The fifth parameter is the 
three-dimensional impact parameter normalized by its error. A ~10% inefficiency in full 
CDC- VXD track reconstruction is partly recovered by also including tracks reconstructed 
in the VXD alone ("VXD-only" tracks). Including VXD-only tracks with at least one hit 
in each of the three VXD layers reduces the tracking inefficiency to ~3% for heavy hadron 
decay products. These tracks are used for vertex charge reconstruction but not for vertex 
momentum or mass reconstruction because of their poor momentum resolution. 

A sample of 53 709 Z° — > bb candidates is then selected by requiring that either hemi- 
sphere in the event contain a seed vertex with mass M > 2 GeV/c 2 , where M is partially 
corrected for missing neutral decay products j^. The fraction of remaining light-flavor events 
is estimated to be 2.9% from the simulation. 
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IV. CHARGE DIPOLE RECONSTRUCTION AND DECAY FLAVOR TAG 



In the last phase of the analysis, we aim to reconstruct both B and D decay vertices 
in a given hemisphere. We can then tag the B° or B° decay flavor based on the charge 
difference between these vertices. The charge difference is expected to be ±2 (0) for decays 
with a charged (neutral) D meson. This charge dipole method relies on the fact that the B 
and cascade D flight directions are very nearly co-linear and the B and D decay points are 
separated along the flight direction. The final vertex reconstruction is carried out separately 
for each event hemisphere and proceeds in two steps. 

In the first step, the 6-hadron flight direction is determined using a set of tracks selected 
by one of two different procedures. The main track selection procedure uses the set of tracks 
(including VXD-only tracks) associated with a seed vertex. An alternative procedure is used 
if no SV is found in the hemisphere (representing 15% of the decays in the final sample). 
In this case, tracks are required to have either >3 VXD hits and p± > 250 MeV/c, or >2 
VXD hits, >23 CDC hits, p ± > 250 MeV/c, a x 2 /DOF < 8 for both CDC and CDC+VXD 
track fits, a first CDC hit with a radius < 39.0 cm, a 2D impact parameter in xy <1 cm, 
and a z distance of closest approach to the IP < 1.5 cm. The latter procedure is useful to 
recover part of the inefficiency of the SV finding algorithm for b hadrons decaying close to 
the IP, a region where the sensitivity to B° s oscillations is highest. The tracks selected with 
either of the two above procedures are then used to estimate the 6-hadron flight direction. 
This is done with a "ghost" track anchored at the IP and initially given an error of 25 /im. 
Each of the selected tracks is individually vertexed with the ghost track and the sum of 
the two-prong vertex fit x 2 values is computed. The direction of the ghost track is varied 
until this x 2 sum is minimized — with the final direction representing the best estimate of 
the 6-hadron line of flight. The ghost track error is then scaled such that the largest single 
two-prong vertex fit x 2 equals 1. 

In the second step, the selected tracks are divided into subsets defining decay vertices 
along the ghost track. This is done using the following iterative procedure. Initially, the 
probabilities for all two-track + ghost track and track + IP candidate vertices are calculated. 
The vertex with the highest fit probability is saved and its track(s) is (are) removed from 
the list of available tracks. If the highest probability vertex is one of the track + IP vertices, 
it defines a new IP vertex. The highest probability vertex then serves as a new candidate 
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vertex for the next step, where again all combinations of the remaining tracks with the new 
candidate vertex, with each other and the ghost track, and with the IP are formed. This 
process continues until all vertex fit probabilities are lower than 1%. A single IP results from 
this procedure. All remaining tracks (not incorporated in the IP or a multi-prong vertex) 
are then combined with the ghost track to form one-prong vertices. One advantage of this 
technique is that it allows for the reconstruction of one-prong decay vertices, which are fairly 
common in both B and D decays. Details about the ability of this technique to reconstruct 
the correct number of vertices can be found in Ref. ^]. 

Hemispheres with exactly two vertices (in addition to the IP) are selected. The (sec- 
ondary) vertex that is closer to the IP is labeled u B n and the (tertiary) vertex farther away 
is labeled U D" . If the D vertex contains two or more tracks, the B vertex is further refined 
by adding a virtual D "track" to those tracks already attached to the B vertex and a new 
vertex fit is performed. The D track is formed using the D vertex position and the net 
momentum direction of the D vertex tracks. The error matrix of the D track is scaled as a 
function of the D vertex mass to partially account for the fact that the D meson is not fully 
reconstructed. 

The vertices must satisfy a series of criteria to guarantee that the B and D vertices 
are well separated and well reconstructed: the B vertex is downstream from the IP and 
has a radius < 2.2 cm, the distance between the B and D vertices (Lbd) satisfies 250 fim 
< Lbd < 1 cm, the D vertex mass is smaller than 2.0 GeV/c 2 (assuming all tracks are 
pions), the ghost track error is smaller than 300 /zm, and the cosine of the angle between 
the nearest jet axis direction and the straight line connecting the IP and the B vertex is 
greater than 0.9. For purposes of flavor tagging (see below), the B and D vertex charges 
are required to be different. Several cuts are applied to enhance the B° s purity and further 
improve the quality of the reconstruction for B° s candidates: the net charge Q to t of all tracks 
associated with the decay chain is required to be zero and neither vertex may contain a 
VXD-only track with p± > 4 GeV/c to make sure that the charge is reliably reconstructed. 
Finally, B% decays are suppressed by vetoing decays including D* + candidates. To this end, 
a veto is applied if slow pion candidates in either B or D vertices have a nD vtx — D vtx 
mass difference less than 0.16 GeV/c 2 . This veto reduces the B® (B^) selection efficiency 
by 5% (13%) and results in a B° s purity increase of 5%. For all data and MC events, we 
remove hemispheres already containing a vertex selected in two complementary analyses, 
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FIG. 1: Distribution of the charge dipole for data (points) and Monte Carlo events (solid his- 
togram). Also shown are the contributions from hadrons containing a b quark (dashed histogram) 
or a b quark (dotted histogram). 

n 

one with fully reconstructed D s |9J and the other with a lepton + D vertex. Such a removal 
is necessary to avoid statistical correlations and facilitate the combination of results from 
different analyses. 

After applying all the above cuts, a sample of 11462 decays remains. The sample com- 
position is 15.7% B + , 57.2% B° d , 15.9% B° s , 10.1% b baryon, and 1.1% udsc, as determined 
from the simulation. The decay flavor is determined by the sign of a vertex charge dipole 
denned as 5Q = Lbd X SIGN(Qd — Qb), where Qb (Qd) is the charge of the B (D) 
vertex. Positive (negative) values of 5Q tag B° (B°) decays. Figure [T] shows the charge 
dipole distribution for the data sample and also indicates the separation between hadrons 
containing b or b quarks in the simulation. The mistag probability for selected B® (B®) 
decays is 22% (24%), as determined from the simulation. The mistag probability depends 
on the charm content of the decay products. For example, the mistag probability is only 
10% for B° s —>■ D S X decays (X represents any state without charm) but 46% for B° s decays 
proceeding through a, b ccs transition, where the latter category accounts for 24% of the 
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selected B° s decays. Detailed studies of B and D vertex charge, multiplicity and separation 
show that the simulation reproduces the data well (see Ref. |5|). 

The proper time is computed using the measured decay length L and momentum p of 
the B° s candidate, T = Lttib/pc. The decay length is defined as the distance between the IP 
and the B vertex. In the particular case of B and D vertices both containing just one track 
(representing 17% of the final sample), the B decay length is reconstructed by averaging 
those decay lengths obtained from the ghost track + B track vertex and the B track + 
D track vertex. The decay length residuals for correctly (incorrectly) tagged B® — > D^X 
decays containing three or more tracks can be described by the sum of a core Gaussian 
with a width of <jl — 78 /im (109 /im) and a tail Gaussian with a width of = 304 fim 
(584 /mi), where the fraction of events in the core Gaussian is 60%. The momentum of 
the B2 candidate is reconstructed with an optimized average of two different methods. The 
starts with tracks associated with the decay and iteratively adds high 
energy (E > 3 GeV) neutral clusters (from the electromagnetic calorimeter) close to the B 
vertex line of flight, until the invariant mass of the decay products is close to the B mass. 
The second method ^j] does not use calorimeter information but relies on the associated 
tracks and the amount of missing transverse momentum to estimate the B momentum. A 
parametrization with the sum of two Gaussians yields a 60% core width of o v jp = 0.07 and 
a 40% tail width of cr p /p = 0.21 for selected B° s decays. 



V. PRODUCTION FLAVOR TAG 



The B flavor at production is determined with a combination of six tags: polarized 
forward-backward asymmetry, jet charge, vertex charge, charge dipole, lepton and kaon 
tags, where all but the first are combined in an "opposite hemisphere charge" tag using 
a series of neural networks. For a detailed description of the polarized forward-backward 
asymmetry Aps see Ref. j^. Briefly, left- (right-) polarized electrons tag b (b) quarks in 
the forward hemisphere, and b (b) quarks in the backward hemisphere. Averaged over the 
acceptance, this yields a mistag probability of 28% for our average electron beam polarization 
of P e = 73%. The opposite hemisphere charge tag employs tracks and topological vertices 
reconstructed in the hemisphere opposite that of the B° s decay candidate. The available 
tags are: (i) the jet charge tag, where tracks are used to form a momentum-weighted track 
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FIG. 2: Distribution of the computed &-quark probability at production for data (points) and 
Monte Carlo events (histograms) showing the b and b components (dashed and dotted histograms, 
respectively). 



charge sum |^|, (ii) the total charge of tracks and charge dipole of a fe-hadron decay, (iii) 
the charge of a kaon identified in the Cherenkov Ring Imaging Detector and associated with 
a 6-hadron decay (if more than one kaon is found, the total kaon charge is used), and (iv) 
the charge of a lepton originating from a 6-hadron decay. These tags are combined using a 
series of neural networks to form an overall production flavor tag characterized by a 6-quark 
probability. The neural networks take into account correlations between the different tags 
as well as vertex charge, mass and decay length dependencies. The corresponding average 
mistag rate is 29%. Finally, the result is analytically combined with the independent b- 
quark probability from the Aps tag- The overall average mistag probability is 22% but the 
information is used on an event-by-event basis, so that a significant fraction of the events 
have a very high production flavor tag purity. The combined tag is 100% efficient. Figure El 
shows the 6-quark probability distributions for data and MC as selected in this analysis. It 
indicates a clear separation between b and b quarks and good agreement between data and 
simulation. 



n 



n 



VI. FIT METHOD AND RESULTS 



Decays are tagged as mixed (unmixed) if the production and decay flavor tags disagree 
(agree). The probability for a decay to be in the mixed sample is expressed as 



e 



-t/r u 



V mixed (t, Am s ) =f u - w p (l - v%) + (1 - w p )w° 

+ ~ ( £ 9dk [(1 - wg)(l + [2w p - 1] cos Am d t) + u&(l - [2w p - 1] cos Am d t) 



Td 



\k=l 



f e~ t/Ts ( 5 r 

+ ^ £ g sk (1 - wf k )(l + [2w p - 1] cos Am s t) + wf k (l - [2w p - 1] cos Am s t) 

* T s \k=l 

g ' t/ TfcjJ- y 

+ fbry [W P (1 ~ W^y) + (1 - W P )w^ y 

T~bry 

+ ^F udsc (t), (1) 

where fj represents the fraction of each fe-hadron type and background (j = u, d, s, bry, 
and udsc correspond to B + , B®, B®, b baryon, and udsc background), F udsc (t) describes 
the proper time distribution of the udsc background, Tj is the lifetime for b hadrons of type 
j, w p is the production flavor mistag probability, and uifo r .y are the decay flavor mistag 
probabilities for B + and b baryons, w^ k and w^ k are the decay flavor mistag probabilities for 
B® and B® with the index k = 1,...,5 representing five different decay final states (D°X, 
D ± X, D S X, charmed baryon X, and final states resulting from a b — ► ccs transition), and 
g d k {.g s k) is the fraction of B® (B®) decays into each of the above final states. A similar 
expression for the probability V U nmixe d to observe a decay tagged as unmixed is obtained by 
replacing the production flavor mistag rate w p by (1 — w p ). 

The quantities in Eq. JI} are determined on an event-by-event basis. The fractions fj 
depend on the angle 9t between the event thrust axis and the beam direction to account 
for the decrease in the Q to t reconstruction purity at high | cos#t|. The mistag probability 
w p depends on cos^, P e , and the opposite hemisphere charge neural network output. The 
mistag probabilities wf k depend on the reconstructed decay length in order to take into 
account the degradation of the charge dipole tag close to the IP. This effect is rather weak 
for Bj, B° s and 6-baryon decays but it is significant for B + decays. Finally, the decay final 
state fractions g dk and g sk are parametrized as a function of M. 

Detector and vertex selection effects are introduced by including a time-dependent effi- 
ciency function e(t) and convolving the above probability density functions with a proper 
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time resolution function TZ(T, t): P m ixed(T, Am s ) = / °° V m ixed(t, Am s ) 7Z(T, t) e{t) dt, where 
t is the "true" time and T is the reconstructed time. A similar expression applies to the 
unmixed probability P unm ixed- Each term of the probability density function is divided by 
a normalization factor given by J 10 ps (P m i X ed + Punmixed) dT, which accounts for the limited 



range of reconstructed proper time considered in this analysis. The reso 
parametrized by the sum of two Gaussian and two Novosibirsk functions 



ution function is 



K(T, t) = 0.36 f Gauss (T, t- a 11 , fi c ) + 0.24 f Gauss (T, t; a u , fi c ) 

+0.24 f Novo (T, t; a 21 , fi t , () + 0.16 f Novo (T, t; a 22 , fi t , (). (2) 

The proper time resolution a 1 ^ is a function of proper time that depends on p, a p , and <Jl, 



2 / j \ 2' 



PC J V P, 



1/2 



(3) 



where the index i = 1 (i — 2) corresponds to the core (tail) component of the decay length 
resolution and the index j = 1 (j = 2) corresponds to the core (tail) component of the 
relative momentum resolution. Slight offsets in the decay length reconstruction are modeled 
by the time-dependent parameters \i c and fi t , whereas an asymmetric tail is modeled by the 
time-dependent parameter (. For each decay, the decay length resolution is computed from 
the vertex fit and IP position measurement errors, with a scale factor determined using the 
MC simulation (scale factors for correctly tagged decays are typically ~1.0 for the core and 
~2.0 for the tail components). The momentum resolution is parametrized as a function of 
the total track energy in each decay, with parameters extracted from the MC simulation. 
The efficiency e(t) is also parametrized using the MC simulation. All parametrizations are 
performed separately for each 6-hadron type and for right and wrong charge dipole tags. As 
a consequence, different resolution functions are used for final states containing one or two 
charm hadrons. 

The study of the time dependence of B° s mixing is carried out using the amplitude 
method |13]. A likelihood fit to the proper time distribution of mixed and unmixed events 
is performed to determine the oscillation amplitude A at fixed values of Am s . That is, in 
the expression for the mixed and unmixed probabilities, one replaces [1 ± cos(Am s t)] with 
[1 ± A cos(Am s t)] and fits for A. This method is similar to Fourier transform analysis and 
has been tested extensively with simulated samples generated at several different values of 
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FIG. 3: Measured B® oscillation amplitude as a function of Am s . The light-grey (dark-grey) band 
is the 90% confidence level allowed region obtained using statistical (total) uncertainties. Values 
of Am s for which the band is below A = 1 are excluded at the one-sided 95% confidence level. 

Am s . The measured amplitude spectrum is shown in Fig. |H1 The measured values are 
consistent with A = for the whole range of Am s up to 25 ps -1 . A signal for B° s mixing, 
corresponding to an amplitude of A = 1, is ruled out for two regions of Am s in the lower 
half of this range, where A + 1.645 a a < 1- 

As a check of the analysis, we performed a fit for the B® oscillation frequency and found 
a value of Am^ = 0.537 ± 0.030(stat) ps -1 , in agreement with the world average value of 
Arrid = 0.489±0.008 ps -1 Q]. Conversely, fitting for the B% oscillation amplitude at a value 
of Am d = 0.489 ps -1 yields A = 0.956 ±0.065, a value consistent with unity as expected for a 
signal. Furthermore, the charge dipole mistag rate is measured directly from the data for B% 
and B + decays and found to agree with MC values within uncertainties (see next section). 
To test the charge dipole tag for B° s decays, we select semileptonic decay candidates tagged 
as mixed and with T < 2 ps to obtain a subsample with a B° s purity of 28%. The charge 
dipole tag is found to agree with the charge of the lepton in (95 ±2)% of the cases in the 
data, to be compared with (94 ± 1)% in the simulation. 
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TABLE I: Measured values of the oscillation amplitude A with a breakdown of the main systematic 
uncertainties for several Am s values. 



Am s 


10 ps- 1 


15 PS" 1 


20 PS" 1 


Measured amplitude A 


-0.517 


0.438 


0.854 


„stat 
a A 


±0.666 


±1.000 


±1.529 


syst 
°A 


+0.217 
-0.269 


+0.444 
-0.284 


+0.254 
-0.364 




-0.123 
+0.162 


-0.149 
+0.136 


-0.194 
+0.170 


B(b — > b baryon) 


+0.028 
n not; 


-0.004 

— U.Uol 


-0.041 

n riQn 
— U.Uoy 


udsc fraction 


+0.002 
+0.005 


-0.023 
—0.002 


-0.055 
—0.024 


Decay length resolution 


+0.031 
-0.041 


+0.003 
-0.002 


-0.022 
+0.028 


Momentum resolution 


+0.042 
-0.126 


+0.024 
-0.183 


+0.056 
+0.003 


Resolution function 


+0.115 
-0.189 


+0.416 
-0.095 


+0.099 
-0.140 


Production flavor tag 


+0.018 
-0.012 


+0.026 
-0.051 


+0.042 
-0.111 


Decay flavor tag 


+0.053 
-0.050 


+0.066 
-0.092 


+0.138 
-0.206 



VII. SYSTEMATICS 

Systematic uncertainties have been computed following Ref. [l3| and are summarized in 
Table H] for several Am s values. The systematic uncertainty takes into account both the 
change in the measured amplitude and in its error 

a syst = A var _ ^nom + ^ _ ^mj A _^ ^ ^ 

a nom 

where A nom (cr^ om ) is the nominal amplitude (error) and A var {<J A ar ) is the amplitude (error) 
obtained for a particular variation of the parameters of the likelihood function. It should be 
noted that the statistical uncertainty dominates for all the Am s values considered. 

Uncertainties in the sample composition are estimated by varying the fraction of udsc 
background by ±20% and the production fractions of B° s and b baryons according to 
0.100 ± 0.012 and 0.099 ± 0.017, respectively [15j. Since the B® purity is one of the most 
important parameters in the analysis, we also varied the and B + branching fractions to 



15 



the five decay final states (D°X, D ± X, D S X, charmed baryon X, and ccsX) by twice the 



uncertainty given in Ref. 



14j . The corresponding branching fractions for B° s decays are var- 



ied by ±20% to account for possible deviations due to SU(3) flavor breaking. Furthermore, 
charge reconstruction uncertainties are evaluated from comparisons between data and MC 
simulation, especially relying on self-calibration tests detailed in Ref. |7J • To account for all 
uncertainties described above, we vary the B® purity by ±13.1%. As a cross-check of the 
fraction of B® mesons in the selected sample, we fit for f s with a x 2 fit to the fraction of 
mixed decays as a function of proper time. A value of 0.159 ± 0.018 is obtained, consis- 
tent with the MC value of 0.159. Other physics modeling uncertainties, which are due to 

n 

uncertainties in the 6-hadron lifetimes and Am^ [15J], are found to be negligible. 

Uncertainties in the modeling of the detector include ±7% and ±10% variations in decay 
length and momentum resolutions, respectively. The decay length resolution uncertainty 
is determined from 3-prong r decay vertices in the data and the momentum resolution is 
confirmed to within 10% by comparing the width of the energy distribution in the data 
with that expected from 6-quark fragmentation . The systematic uncertainty due to the 
particular choice of proper time resolution function is evaluated by varying core fractions 
and offset corrections. 

The production flavor tag systematic error accounts for the uncertainties in the forward- 
backward asymmetry and opposite hemisphere charge tags. Dominant uncertainties for 
the forward-backward asymmetry tag are the electron beam polarization P e = 0.73 ± 0.02 
and the parity-violating parameter = 0.935 ± 0.040. The uncertainty in the opposite 
hemisphere charge tag is extracted from events in which both hemispheres have such a tag. 
These double-tagged events provide a direct measurement of the average mistag rate, which 
agrees very well with the MC value and has an uncertainty of 0.008. The total uncertainty 
in the average production flavor mistag probability is estimated to be 0.008. Furthermore, 
we allow for a deviation in the shape of the 6-quark probability distribution for the data. 
The corresponding uncertainty is dominated by a 7% uncertainty in the slope of the 6-quark 
probability as a function of opposite hemisphere charge neural network output, as extracted 
from double-tagged events in the data. 

The decay flavor mistag probability is varied by ±7% for all b hadrons. This uncertainty is 
determined from \ 2 fits to the B® and B + mistag probabilities in the measured dipole-tagged 
forward-backward asymmetry and the fraction of mixed decays as a function of proper time. 
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These fits confirm that the mistag rate for B® and B + decays is correctly modeled in the 
simulation. To account for a possible deviation from SU(3) flavor symmetry, an additional 
7% uncertainty is assigned to decays. This value is obtained assuming an additional 
20% uncertainty in the branching fraction for B° s decays proceeding via a b — > ccs transition 
(since these dominate the mistag rate). 

Including the systematic uncertainties, the following ranges of B® oscillation frequencies 
are excluded at the 95% C.L.: Am s < 4.9 ps _1 and 7.9 < Am s < 10.3 ps _1 . That is, the 
condition A + 1.645 a a < 1 is satisfied for those Am s values. The sensitivity to set a 95% 
C.L. lower limit, defined as the Am s value below which 1.645 a a < 1, is found to be 8.7 



VIII. CONCLUSIONS 

We have performed a search for B° s -B® oscillations using a sample of 400 000 hadronic 
Z° decays collected by the SLD experiment. The cascade structure of 6-hadron decays is 
reconstructed and the B flavor at decay is tagged with a novel charge dipole technique 
that relies heavily on the excellent resolution of the CCD pixel vertex detector. The B 
flavor at production is tagged using the polarized electron beam as well as track and vertex 
information from the hemisphere opposite that of the B® candidate. With a final sample of 
11462 decays, we achieve a sensitivity of 8.7 ps _1 and exclude the following values of the 
Bg oscillation frequency: Am s < 4.9 ps -1 and 7.9 < Am s < 10.3 ps _1 at the 95% C.L. 

Combining these results with those obtained from a sample of fully reconstructed D s 
mesons jj] yields a combined SLD sensitivity of 11.1 ps -1 and excludes B° s oscillation fre- 
quencies < 10.7 ps -1 at the 95% C.L. This result confirms previous studies performed by 
the ALEPH, CDF, DELPHI, and OPAL Collaborations di]], and is among one of the most 
sensitive to date. 
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